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Abstract. In this work I investigate the statistical properties of a huge catalog of 
closely interacting pairs of dark matter haloes, extracted from the Millennium Simu- 
lation database. Only haloes that reach a minimum mass > 8.6 x 10 lQ M G h~ l (corre- 
sponding to 100 particles) are considered. Close pairs are selected if they come within a 
critical distance d a it- I explore the effects of replacing d a it = 1 Mpc/z -1 — > 200kpc/i~' 
on the evolution of separations, lifetimes, total masses and mass ratios of these pairs. 



1. Introduction 

Mergers of galaxies play a fundamental role in essentially all modern theories of galaxy 
formation. They are believed to determine the morphology of a galaxy, drive its star 
formation and even activate its nuclear supermassive black hole. This has motivated 
the development of very detailed numerical simulations of merging galaxies (often in- 
volving only two galax ies in isolation). Unfortunately, with very few exceptions (e.g., 



Tonnesen & Cenll201 1 ). these simulations typically concentrate primarily on the post 
merger aftermath, often neglecting the early stages of interaction. 

Observations, on the other hand, tell us a different story. Surveys like the SDSS 
and zCOSMOS have confirmed, in large numbers, that galaxies in pairs te nd to be bluer, 
have t heir star formation enh anced, and are more likely to be active (e.g. jEllison et al. 



201 ll : IsTiverman et al.l 1201 lT) . Moreover, the discovery rate o f binary quasa rs has ac 



celerated to unprecedented levels in the last few years (e.g.. iLiu et al.ll2011l) . For this 



reason, it is vital to re-focus our attention to the early stages of galactic interactions. 



2. The Halo-Pair Catalog 

This work is based on the results of Moreno (2011), where I use the publicly avail 



able Millennium Simulation database (ISpringel et all 120051) to construct a very large 



catalog of closely-interacting pairs of haloes. First, pairs are selected if they come 
within a critical distance cf cr it. Also, I require that all haloes reach a mass of at least 
> 8.6 x 10 10 M G h (corresponding to 100 particles). Below I explore the effects of re- 
placing d CT [ t = 1 Mpc/j -1 — > 200 kpc h (these two sets are named the 'Close Set' and 
the 'Very Close Set' respectively). The true physical separations at which interaction- 
induced phenomena are instigated may depend on the environment and on the actual 
phenomenon we care about. For this reason, at this preliminary stage, I choose the 
simplest possible criterion. More sophisticated conditions can always be implemented, 
and the corresponding catalogs can then be extracted from this primary set. 
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Fi gure 1. Pairs of branches extracted from the merger tree presented in Figure 1 
of Mo renol (1201 lh . Red solid circles represent mergers, while dotted parallelograms 
refer to close interactions. LEFT-PANEL: Branches B and C merge (love). RIGHT- 
PANEL: Branches C and D interact for a while, but split eventually (deception). 



To describe the evolution of each pair, I propose a scenario based on five-stages: 
initial, entry, closest, final and fate. The initial, closest and final stages are those at 
which the two haloes are identified in the simulation for the first time, at their closest 
distance, and for the last time respectively. The entry stage is the first time the pair is 
at a distance lower than <i cr i t . The fate stage comes in two flavors: either the two haloes 
merge or not. The latter situation happens because either they do not have enough time 
to merge by the present time or one of the haloes is absorbed by a third halo. This 
process is analogous to a situation that occurs commonly in real life. Suppose you meet 
a person and the two of you date for a while. The following two things can happen: 
either you marry this person eventually (love) or someone else takes this person away 
from you (deception). In the simulation, the culprit of this deception is usually a very 
massive halo in the vicinity. Figure 1 illustrates these two situations. 



3. Results 

For the sake of readability, I present all my results together in Figure 2. The corre- 
sponding colors and line-styles are indicated in the key. To show the effect of shrinking 
J cr it = 1 Mpc h — > 200 kpc h~ , I split each of the four panels in two: the upper por- 
tions refer to the Close Set, while the lower portions correspond to the Very Close Set. 
Before going into details, the first thing to notice is that when a smaller d CI i t is reduced, 
halo pairs that do not merge are less likely to be selected (from 53% — > 42%). 

Halo pairs usually begin at very high redshift with separations ~ 10 Mpc h~ x . 
They cross the critical threshold, reach their minimum separation, and are identified as 
two distinct objects for the final time. While a fraction of these pairs end up with sub- 
critical separations, many of these have super-critical final separations. This fraction 
increases as <i cr i t is reduced (from 45% — > 57%). In other words, using a smaller 
critical distance selects pairs that end up with more elongated orbits. For those haloes 
that never merge, shrinking d cr n actually augments the fraction of super-critical fate 
separations (from 40% — > 80%). In other words, while some super-critical pairs at the 
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Figure 2. Results. Each of the four panels is split in two: the upper portions refer 
to d CS ;it = 1 Mpc h~ l (Close Set), while the lower portions refer to af cl -i t = 200 kpc h~ l 
(Very Close Set). UPPER-LEFT: Separation. UPPER-RIGHT: Lifetime. LOWER- 
LEFT: Total mass m\ +rri2. LOWER-RIGHT: Mass ratio p = m\/m2 (with m\ < mi). 



final stage come back and merge, others separate permanently. Using a smaller critical 
distance enhances the contribution of cases where this splitting is more violent. 

The total lifetime of the pair is the sum of the approaching time (from initial to 
entry) and the interacting time (from the first to the last time the separation is below 
cfci-it). The total duration of the pair can be very long, almost always larger than a few 
Gigayears. In fact, most of this time is spent in the approaching phase. The distribution 
of the total and infall times is largely insensitive to the choice of critical distance. How- 
ever, the interaction times tend to be shorter when d crit is reduced (e.g., the fraction of 
pairs with interactions shorter than 1 Gyr goes from 40% — > 60%). 

Haloes in pairs typically begin with masses just above the resolution of the sim- 
ulation. As time progresses, the sum of the masses tends to increase at scales > 
3 x 10 13 M G /j _1 , otherwise they decrease. In other words, this particular mass scale 
marks the transition at which the evolution is dominated by tidal shredding as the pair 
sinks within a larger host. The great majority of pairs remain in the small-mass regime 
throughout their existence. Reducing c? cr it increases the fraction of pairs in the massive 
regime (from 3% — > 8% during the interaction phase). In contrast, the fate mass dis- 
tribution is very different from the others because it is largely determined by external 
haloes in the vicinity. Notice that the fate total masses tend to be extremely large! 

Since the initial masses are usually near the resolution of the simulation, pairs tend 
to be born with similar masses (and thus with mass ratios p < 1; e.g., in the 'major' 
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pair regime). Throughout their evolution, pairs tend to stay more or less major; and 
only at the interaction stages do the more discrepant 'minor' pairs begin to appear. The 
importance of minor pairs becomes more prevalent as d C n t decreases (from 19% — > 33% 
for fi < 1/10). Nevertheless, minor pairs never dominate the bulk of the population. In 
contrast, notice that fate mass ratios tend to be extremely small (i.e., very discrepant)! 

In particular, the total-mass and mass-ratio distributions point to the fact that at 
small scales (probed with the smaller d CT i t ), the fraction of pairs involving a very mas- 
sive member increases. Moreover, these results in conjunction also appear to unveil 
the real culprits responsible for splitting close pairs: very massive haloes in the vicinity 
(usually the central halo in a group or some other massive secondary halo) absorb one 
of the members of the pair before these two have a chance to merge. And as smaller 
scales are probed, the effects of these vicious culprits become even more evident. 

4. Conclusions 

Using the Millennium Simulation, I analyze the behavior of a very large set of halo 
pairs selected by a proximity criterion. As this critical distance is reduced, (1) more 
mergers are pre-selected, (2) pairs end up with more elongated orbits, (3) interactions 
tend to be shorter, and (4) non-mergers split more violently. Using a smaller <i cr it probes 
smaller scales, enhancing the importance of pairs involving a massive member and the 
prevalence of massive neighbors capable of splitting pairs. 

With these results, I wish to highlight the following important lessons: 

• Extra care must be taken when using close galaxy pairs as proxies for merg- 
ers. This subtlety is important even in physical three-dimensional space! 

• Systems with two merging galaxies in isolation are just an approximation. 
In reality, the universe can nurture these duos, or split them altogether! 

Next I will use this catalog to explore the symbiosis between interacting galaxies 
and binary quasars. Ultimately, this will set up realistic initial conditions for black hole 
mergers. This work is just a puzzle piece of an ambitious long-term research program 
centered on the evolution of galaxies, supermassive black holes, and their environment. 
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